Gravitational Radiation Generated by Cosmological Phase Transition Magnetic Fields 
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We study gravitational waves generated by tiie cosmological magnetic fields induced via bubble 
collisions during the electroweak (EW) and QCD phase transitions. The magnetic field generation 
mechanisms considered here are based on the use of the fundamental EW minimal supersymmetric 
(MSSM) and QCD Lagrangians. The gravitational waves spectrum is computed using a magne- 
tohydrodynamic (MHD) turbulence model. We find that gravitational wave spectrum amplitude 
generated by the EW phase transition peaks at frequency approximately 1-2 mHz, and is of the or- 



der of 10" 



10 ; thus this signal is possibly detectable by Laser Interferometer Space Antenna 



(LISA). The gravitational waves generated during the QCD phase transition, however, are outside 
the LISA sensitivity bands. 

PACS numbers: 98.70.Vc, 98.80.-k,98.80.Cq 
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I. INTRODUCTION 

There are various mechanisms that might generate 
gravitational waves in the early Universe. For reviews 
see P, 0, 3. A well-known one is the generation of 
gravitational waves during the cosmological EW or QCD 
phase transitions [1, S ll, 0|- These mechanisms in- 
clude bubble wall motions and collisions if the phase 
transition is first order [§| , as well as cosmological mag- 
netic fields and hydrodynamical or MHD turbulence 

A brief overview of our present work is: 

• We use basic EW and QCD Lagrangians to derive 
the magnetic fields created in the plasma during 
the EWPT and the QCD phase transitions. 

• In the framework of the standard MHD theory the 
magnetic field produces turbulence in the plasma. 
The value of the Alfven velocity, va, is a key pa- 
rameter when considering the generation of gravi- 
tational waves. 

• From the parameters found via our theory of the 
EW and QCD phase transitions, we use the formal- 
ism of the gravitational waves generation by MHD 
turbulence to estimate the peak frequency and am- 
plitude of the gravitational waves produced during 
these cosmological phase transitions.. 

The direct detection of the relic gravitational waves 
will open the new prospects to understand the physical 
processes in the early Universe [ij. The main objective 
of the present paper is to study if the gravitational waves 
produced by the magnetic fields created during cosmolog- 
ical phase transitions could be detected by the current 
and/or nearest future missions. To be observable the 
gravitational waves signal must satisfy two conditions: it 



must be within the observation frequency bands, and its 
amplitude should exceed substantially the instrumental 
noise (for the stochastic backgrounds signal to noise ratio 
(SNR) must be taken to be 5) Our present study is 
to determine if the gravitational waves produced either 
by the EW or the QCD phase transitions might be de- 
tectable by LISA, whose sensitivity reaches maximum at 
frequencies 1 — 100 mHz Q. 

In order to produce a detectable gravitational wave 
signal the cosmological phase transition must be first or- 
der, so that bubbles of the new vacuum nucleate within 
the false vacuum |^ at a critical temperature. Other- 
wise, there is a crossover transition. It has been shown 
that with the Standard EW model there is no first or- 
der phase transitions [T6| . and there is no explanation 
of baryogenesis. However, there has been a great deal 
of activity in the MSSM extension of the standard EW 
model [131 • In this case a MSSM having a Stop with a 
mass similar to the Higgs mass leads to the first order 
EW phase transition. This EW MSSM theory is consis- 
tent with baryogenesis [l^l- On the other hand, recent 
lattice QCD calculations have shown flO^ that the QCD 
phase transition is a first order cosmological phase tran- 
sition, with bubble nucleation and collisions. 

The EW phase transition is particularly interesting for 
exploring possible cosmological magnetic fields since the 
electromagnetic field along with the and Z fields are 
the gauge fields of the Standard model. For the QCD 
phase transitions the electromagnetic field is included 
in the Lagrangian through coupling to quarks. In both 
cases the magnetic field with large enough energy density 
can have different cosmological signatures. In particular, 
the big bang nucleosynthesis (BBN) bound on the mag- 
netic field energy density is: pB < O.lprad [13]) gi'^ing 
= Ps/pcr < 2.4 X 10~^/io^, where flB is the energy 
density parameter at the present time, pcr is the present 
critical energy density (i.e. the total energy density for 
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the flat Universe), and ho is the current Hubble parame- 
ter in the units of lOOkm/sec/Mpc. This corresponds to a 
hmit on the effective magnetic field comoving amplitude 
~ ^SirpB < 8 X 10^'' Gauss. ^ Similar Hmits can be 
obtained through the available data of Laser Interferom- 
eter Gravitational Observatory (LIGO), if it is assumed 
that the primordial magnetic field generates relic gravita- 
tional waves via its anisotropic stress [2l| . Stronger lim- 
its on the primordial magnetic field are provided by the 
cosmic microwave background (CMB) anisotropy data, 
B'^^ < 5 X 10~^ Gauss, see Ref. [l^ for a review and 
references therein. 

Several studies have been performed to estimate the 
gravitational wave signal from the first order EW phase 
transition in MSSM and nMSSM, see Refs. Using a 
model similar to the tunneling from a false to a true vac- 
uum [1] it was argued that gravitational waves generated 
by the EW phase transition are possibly detectable 0. 
However, based on a similar model, it was concluded that 
the magnitude of the gravitational waves generated in a 
MSSM model of the EW phase transition have a lower 
amplitude than that required to be detected by LISA, 
while in a nMSSM model the signal is larger 24 1. 

In the present work we re-address the gravitational 
waves generated by the primordial magnetic fields cre- 
ated during the cosmological phase transitions. There are 
several models of phase transitions resulting of the mag- 
netic field production ^25j . For reviews and references, 
see Refs.[26l.[27|. Our calculations for the EW phase tran- 
sition are based on previous studies of the magnetic field 
generation by nucleation [28, 29] and collisions [30, 31] of 
the EW bubbles. In these studies the basic MSSM EW 
Lagrangian is used. Similarly, to generate magnetic fields 
during the QCD bubble collisions [H, the basic QCD 
Lagrangian is used [33]. For the QCD phase transitions 
the bubble walls are composed mainly of the gluonic field. 
After the collision of two bubbles interior gluonic wall is 
formed, resulting in a magnetic wall production (due to 
coupling of quarks within nucleons to the gluonic wall 
causing alignment of nucleon magnetic dipole moments) . 

The magnetic fields generated by the cosmological 
phase transitions lead to MHD turbulence. In the present 
work we make use of the basic MHD formulation [S^] and 
determine the main parameter for MHD turbulence - the 
Alfven velocity, va, for the magnetic fields produced both 
by the EW and QCD phase transitions. The MHD tur- 
bulence model described in details in Refs. [l^ is 
then used to calculate the gravitational waves produced 
by MHD turbulence present during EW or QCD phase 
transitions. In contrast to the previous works P.[l3l[36j. 
we do not parameterize the gravitational wave signal in 
terms of certain phase transitions parameters, but use 
our solutions based on the fundamental EW and QCD 



Lagrangians. On the other hand, our results are model 
dependent in the sense that they use the value of the bub- 
ble wall velocity Vb = 1/2, found in Ref. [l^. However, 
the value of Wf, = 1/2 was derived using the fundamental 
EW MSSM theory for bubble nucleation rather than a 
model (see below). 

The outline of our paper is as follows. In Sec. II we re- 
view the main MHD equations, the origin of turbulence, 
and discuss the main assumptions used in our present 
work. In Sec. HI we discuss the EW and QCD equations 
used for deriving the magnetic fields produced by the two 
cosmological phase transitions, and give our results. In 
Sec. IV we study the gravitational waves generated by 
these phase transitions driven by MHD turbulence. In 
Sec. V we present and discuss our work. In Sec. VI we 
give our conclusions. 



II. MHD TURBULENCE MODEL 

Our objective is to derive the gravitational wave energy 
density, pcw produced by the magnetic field, B, created 
during EW and QCD phase transitions. In the first sub- 
section section we discuss the basic MHD theory leading 
to turbulence, and in the second subsection review the 
MHD turbulence model used for our cosmological appli- 
cations. 



A. MHD Turbulence and the Alfven Velocity 

An essential aspect of our model is that the magnetic 
field created after bubble collisions couples to the plasma 
and creates MHD turbulence, which then generates via 
the anisotropic stress the gravitational waves. In other 
words, the magnetic energy density pB is transformed to 
the gravitational waves energy density paw- Since the 
coupling between the magnetic and gravitational waves 
energy occurs through the Newton gravitational constant 
G, due to the small value of G the efficiency of the gravi- 
tational wave production directly from the magnetic field 
is too small, see also Refs. [l^, but even accounting for 
the small efficiency, the gravitational wave signal can be 
possibly detectable. 

To show coupling between an initial phase transition 
generated magnetic field to the plasma we give here the 
basic MHD equations for an incompressible, conducting 
fluid Isi 



d_ 

drj 



+ (v • V) - vV' 



(b • V)b - Vp -f f , (1) 



d ^2 
ari 



^ In what follows we use the natural units c = ?i.= fes = l and 
MKS system for the electromagnetic quantities. 



b = -(v V)b + (b- V)v,(2) 
Vv = Vb = , (3) 

where "q is the conformal time, v(x, rf) is the fluid veloc- 
ity, b(x, ?/) = B(x, 7y)/v47rw is the normalized magnetic 
field, f(x, ry) is an external force driving the flow, v is 
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the comoving viscosity of the fluid, r^ro is the comoving 
resistivity, and w= p + p, p, and p the enthalpy, energy 
density, and pressure of the plasma. 

The coupling of the magnetic field to the plasma (see 
Eqs. (HI)-® leads to Alfven turbulence development with 
the characteristic velocity, va [13, HI] , 

B iSpB 
■"A = -7= = a/ ^ • (4) 

Since both pB and p — piad scale as l/a*(i), (where 
a{t) is the cosmic scale factor), if there is no damping 
of the magnetic field (additional pB temporal depen- 
dence), the value of va is not affected by the expan- 
sion of the Universe, va is an essential parameter in 
the generation of gravitational waves in the magnetized 
turbulent model considered here (for details of the hydro- 
turbulence model, see Ref. [l^l)- To develop the turbu- 
lence picture it is assumed apriori that the EW or QCD 
bubble collisions lead to the vorticity fluctuations, i.e. 
the presence of the kinetic turbulence, with a character- 
istic velocity (associated with the largest size bubble) vq. 
Equipartition between the magnetic and kinetic energy 
densities implies va — vq.'^ We also define the energy 
density of the plasma to be equal of the radiation en- 
ergy density, prad, and prad at the moment of the phase 
transition with temperature T^, is given by 

Prad(at T,) = — g^T,)" , (6) 

where gi, is the number of relativistic degrees of freedom 
at temperature T^, [39,]. Using the BBN bound on the 
total magnetic energy density pB < O.lprad, the Alfven 
velocity must satisfy va < 0.4. 

The Afven velocity va as well as the bubble kinetic 
motion velocity vq can be related to the phase transi- 
tions parameters, apT, the ratio between the latent heat 
and the thermal energy, and the efficiency, kpt, which 
determines what part of the vacuum energy is trans- 
ferred to the kinetic energy of the bubble motions as 
opposed to the thermal energy. Ref. Q presents the 
estimate for the largest size turbulent bubble velocity, 
VA — Vq — ^ kptQ;pt/(4/3 + kptQ!pt), and thus ac- 
counting for va — Vq^ BBN bound leads kpt^pt < 0.2. 



Accounting for the stochastic nature of the magnetic field, the 
Alfven velocity is scale dependent, and if it is not specified, va 
is associated with the largest size magnetic eddy. On the other 
hand, va can be expressed in terms of the magnetic field comov- 
ing amplitude as, 

.,^4xl0-*f^,^ \{^y'^\ (5) 

VlO-»Gauss ) VlOO/ 

Here we used Eqs. Q and The MHD turbulence descrip- 
tion presumes that the magnetic turbulent energy density is sat- 
urated when the Alfven velocity reaches the kinetic velocity of 
the plasma, i.e. va — "fo- 



B. Direct Cascade and Magnetic Energy Density 

After being coupled to the fluid, the magnetic field en- 
ergy density injected into the plasma at characteristic 
comoving length scale Aq (which is, of course, inside the 
comoving Hubble radius at the moment of the field 
generation). The magnetic field energy density due to 
coupling with the fluid motions is re-distributed spatially 
through the following regimes [11], fc < /cq (^o = 27r/Ao), 
ko < k < ko (with kjy — 2'k/Xbi the wavenumber cor- 
responding to the magnetic field damping scale due to 
the plasma viscosity) and k > ko- Inside so-called in- 
ertial range, kg < k < ko, the selective turbulence de- 
cay occurs and the magnetic energy fiows from the large 
to the small scales according to the direct cascade Kol- 
mogoroff law j40| . resulting in the magnetic field spec- 
tral energy density EM{k) oc fc^^/"^, while at large scales 
[k < ko) the free turbulence decay takes place, lead- 
ing to the magnetic field spectrum EM{k) oc /c"^, with 
ttT > 4 [4l|. The initial Batchelor spectrum Q!t = 4 
can be transformed via the non-linear processes to the 
Kazantzev spectrum, ax = 3/2 ]42i] . Another numeri- 
cal realization of large scale turbulence can be the white 
noise spectrum, ut = 2, (Saffman) spectrum [i^. Us- 
ing the fact that the total energy density of the magnetic 
field Em = Jq dk EM{k) can not be larger than the ini- 
tial magnetic energy density, pB, it is straightforward to 
obtain the maximal allowed values for EM{k), and get 
the magnetic field limits at large scales [13 ]. 

In our present work MHD turbulence is created by bub- 
ble collisions. It is physically justified to assume that the 
typical injection scale of the magnetic energy is associ- 
ated with the phase transition largest bubble sizes, which 
are given by the bubble wall velocity vt and phase tran- 
sition time-scales f3~^ for the EW and QCD phase tran- 
sitions (determined through the bubble nucleation rate) , 
i.e. 7 = Xo/Xh = VbiP/H^,)-'^ < 1 fl. With these as- 
sumptions one finds for the Kolmogoroff power law and 
the wave numbers [l3| 

EM{k,t)^CMS^^^k-^^^ (7) 

over the range of wave numbers ko < k < ku, where k]j = 
koRe^^^, Cm — 1, -R ^ 1 is the turbulence Reynolds 
number at the temperature T^, and e = (2/3) •^/^/cqv^ is 
the comoving magnetic energy dissipation rate per unit 
enthalpy. 

We consider only the inertial (direct cascade) range 
due to following reasons [11]: i) helicity vanishes (for 
the effects of initial kinetic or magnetic helicity, see 
Refs. [12] •) in our symmetric treatment of bubble col- 
lisions for the phase transitions; ii) we presume that the 
contribution from the large scales {k < kg) into the gravi- 
tational waves signal will not exceed significantly (or even 
being smaller) than that which comes from the inertial 
range, due to the free decay of non-helical turbulence 
and the magnetic energy density small amount presence 
at large scales (we will address this issue in the separate 
work [45|). Another important assumption that we make 
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^ within the inertial range for the 
= fco, as a consequence, the smallest 



is equipartition between the kinetic and magnetic energy 
densities, which allows us to use the direct analogy be- 
tween the hydro and magnetized turbulence. 

The direct cascade turbulence is characterized not only 
by the spatial structure (fc-dependence), but it is impor- 
tant to take into account the time dependence of the 
turbulent quantity correlations. First, we assume that 
the source lasts enough time to allow us to consider the 
developed turbulence, so we do not include in our calcu- 
lations the pulse-Hke source [l^. To compute the direct 
cascade duration time we use the fluctuation time decor- 
relation function ri{k) ~ y/2Tr{k/ko)'^^^TQ^ [i^, where 
To = lo/vo is the largest turbulent eddy turn-over time. 
Note lo = vi,P~^ is the physical length scale of the largest 
size eddy (bubble). As a result the time dependence of 
the magnetic field two-point correlation function within 
the inertial range is given by the function f{ri[x),T) 
[4^ . The function f{r]{x),T) is such that it becomes 
negligibly small for r ^ l/r](k), and from the dimen- 
sional analysis is adopted /[?7(fc), t] — exp[— 7r77^(fc)r^/4] 
[46l |. It is clear that the temporal decorrelation func- 
tion r]{k) oc (fc/fco)^^^ reaches its maximum equal to 

maximal size eddy, k = 
eddies are decorrelated first, and the turbulence cascade 
time is determined by the largest size eddy turn-over time 
tq . As a result the comoving (measured today) peak fre- 
quency of the induced gravitational waves is given by 
/peak ^ va/Xo (see below). 

The proper consideration of the temporal decorrelation 
leads to the fast damping of the gravitational wave sig- 
nal amplitude for frequencies / ^ /peak, i-e. larger than 
that associated with the direct cascade turbulence in- 
duced peak frequency, [ll,[l3l- Several previous studies, 
Refs. |Tol . [Tl| . did not account for the temporal expo- 
nential decorrelation function, and as a result the shape 
of the gravitational wave at high frequencies was given 
by power law, without having steep exponential damp- 
ing ^is'l . We also underline that our description does not 
apply for any pulse-like sources [15i] . In the former case, 
the characteristic comoving gravitational wave peak fre- 
quency is determined by I/Aq [Ti| . so it is higher than 
that in the case considered here. 

Another consequence of the temporal decorrelation is 
that the turbulence cascade time-scale is much shorter 
than the Hubble time scale, and thus we are allowed 
neglect the expansion of the Universe during the grav- 
itational waves generation process [i^. We account 
for the expansion of the Universe only when computing 
the gravitational wave amplitude hc{f) measured today 
(or the corresponding spectral energy density parameter 
51gw(/)), as a function of the linear frequency / mea- 
sured today. 



III. MAGNETIC FIELDS GENERATED BY EW 
AND QCD PHASE TRANSITIONS 

A. Electroweak Phase Transition 

First we review the magnetic field created during EW 
phase transition. In a suitable MSSM model the EW 
phase transition is first order, which results in bubble nu- 
cleation and coUisions. The MSSM EW theory for bubble 
nucleation is developed in Ref. [28i], a Weinberg-Salam 
model with all supersymmetric partners integrated out 
except the Stop, the partner to top quark, has the form 



MSSM 



£1 + £2 + £3 

+leptonic and quark interactions (10) 



^r-W,~'-B,m^-V{^) 
9s 



with 



d,,W:-d.W;-9e,,kW^W!: (11) 



B^v = d,_iBy — d^B^ , 



where the W\ with i = (1,2), are the W+,W~ fields, 
is an SU(3) gauge field, ($, $5) are the (Higgs, right- 
handed Stop fields), (t^X") are the (SU(2), SU(3)) gen- 
erators, and the EM and Z fields are defined as 



A em 



1 



+ g" 



ig'W;i+gB,) 



\/9 + 9 



(12) 



The parameters used here are g — e/smOyy — 0.646 
and g' = g tan^vy = 0.343. The equations of motion 
were solved using an SU(2) I-spin ansatz for the gauge 
fields, and it was found that the nucleation velocity of 



related through 



hc(f) = 1.26 X io-« (^^^^) kcgwc/)]'^', 



(8) 



where / is the linear frequency measured today, and it is given by 
/ = {2it)~^u!, with w = (a*/ao)tiJ*, where oj* is the gravitational 
wave frequency at the time of generation. Due to the Universe 
expansion the freely propagating gravitational wave amplitude 
and frequency are rescalled by a factor 



X 10" 



ao 



100 GcV \ 



100 \ 

g* J 



1/3 



(9) 



The amplitude and the energy density of the gravitational wave 



where a* and ao are the scale factors at the time of generation 
and today, respectively. See Ref. [l|] for the definitions and dis- 
cussions on gravitational wave direct detection experiments. 
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the bubble walls, Vb = 1/2. This is a very important pa- 
rameter for our present work. Using EW MSSM theory, 
we find Vb to be larger than found using models for EW 
bubble walls [i?!. 



Starting from the Lagrangian given by Eq. (jlOp for 
the MSSM EW phase transition, from which bubbles are 
formed, to compute bubble collisions we assume that the 
Stop is integrated out, so the term is not included. 
The magnetic field created during the collision of two 
bubbles was estimated by Refs. [sO, [lH. Our model is 
based on the assumption that the final magnetic field 
is created by the final two bubbles colliding. The time 
scale is from 10~^^ to 10"^'' seconds, and the criti- 
cal temperature (EW phase transition energy scale) is 
~ Mh, with the Higgs mass Mh ^ 110 to 130 GeV 
ps'l . The magnitude of the magnetic field created in this 
final collision of two bubbles, B^^^\ is found in our new 
calculations based on Ref. [3l|, but with the coUid- 
ing bubbles having larger overlap. That is the bubble 
collision has been followed to the stage of overlap of the 
bubbles such that the B-field in the Universe just after 
the EW phase transition has been determined. The re- 
sult for the B-field is 



^ lOM^ 
= 6.4 X lO'^M.^ 



W,80 



GeV^ . (13) 



where Mw,80 — Mw/SOGeV is normahzed W-boson 
mass. Even though the above estimate is given for two 
bubbles collision, we extend our consideration presuming 
that there is a continuous creation of the magnetic field 
through the bubble collisions, and the total magnetic en- 
ergy density released during the MSSM EW phase tran- 
sition can be approximated by — (B^^^^)"^ /8tt. In 
reality of course there are many bubble collisions,'* which 
leads to the establishment of the stochastic nature of the 
resulting magnetic field. The coupling between this ini- 
tial field with the fluid results in MHD turbulence [IJ- 
Using Eqs. ([4])- ([6]) we find that the Alfvcn velocity is 



(EW) 



0.27M?, 



W.80 



VlOOGeVy VlOO/ 



-1/2 



(14) 



B. QCD Phase Transition 

Recent lattice QCD studies have shown that the QCD 
phase transition is first order [l^. In Ref. [13] the mag- 
netic field created by the QCD phase transition bubble 
collisions was derived. A gluonic wall is created as two 
bubbles collide, and a magnetic wall is formed by the 
interaction of the nucleons with the gluonic wall, with 
electromagnetic interaction Lagrangian 



(15) 



where ^ is the nucleon field operator and A^™ is the elec- 
tromagnetic 4-potential. In Ref. it was shown that 
the interaction of the quarks in the nucleons with the glu- 
onic wall alignes the nucleons magnetic dipole moments, 
producing a B-field orthogonal to the gluonic wall. 

Using an instanton model, for the gluonic instanton 
wall oriented in the x-y direction one obtains for = 
^(QCD) ^jjg wall, of thickness ( 



CAqcd 2M„ 



(16) 



where Aqcd is the QCD momentum scale. A similar form 
had been derived earlier using a domain wall model [sol ] . 
The value for B'^^'^^^ was found to be 



(QCD) 



0.39-A|cD - 1-5 X lO^^GeV^ . (17) 



Eqs. (HI)-® with g,=15, T*=0.15 GeV give 



V 



(QCD) 



8.4 X 10"^ 



(18) 



IV. 



GRAVITATIONAL WAVES GENERATED 
BY MAGNETIC FIELDS 



In this section we derive the strain amplitude mea- 
sured today, he (/) , of the gravitational waves generated 
by MHD turbulence developed during the EW or QCD 
phase transitions. The gravitational waves are generated 
by the transverse and traceless part, Sij, of the stress- 
energy tensor, Tij [sij . 



1 



5,,(x,t) = T,,(x,<) - -%r^-(x,t) . 



(19) 



For the EW phase transition with ~ 100 and ~ 100 

GeV, one finds that w^^'' ~ 0.27, which is below the 
BBN bound discussed above. On the other hand the 
value of VA — 0.27 means that 7.4% of the radiation en- 
ergy density is in the form of the magnetic energy density. 



* It is clear that within the Hubble radius we have several areas 
where the magnetic field is generated. The size of the colliding 
bubble determines the correlation length of the magnetic field, 
as we noted above. 



When considering gravitational waves generated by the 
magnetic field, the source Sij is associated with the mag- 
netic anisotropic stress Q. On the other hand, if the 
duration of the source is short enough when comparing 
with the Hubble time at the moment of the generation, 
H^^, we can neglect the expansion of the Universe and 
the gravitational wave generation is described by the sim- 
plified equation and solution 

V'/iy(x,i)- — /iy(x,t) = -16^G5.y(x,i) 



(20) 
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Here hij (x, t) is the metric tensor perturbation which sat- 
isfies the following conditions: ha — and dhij/dx^ — 0. 

To derive the energy density of the induced grav- 
itational waves, accounting for the stochastic nature 
of the magnetic turbulence source, and as a conse- 
quence the stochastic nature of the gravitational sig- 
nal, we must compute the autocorrelation function 
{dthij (x, t)dthij (x, t + t)) /327rG which can be expressed 
through the two-point correlation function of the source 
Rn^J{Lr) = (5y(x',t)5„(x",< + t)), where | = x" - x' 



[Til . As the calculations performed in Refs. [H, [T3| 
show the gravitational energy density, pcw at the mo- 
ment of generation is given by the duration time tt and 
the Fourier transform of Rijij tensor, Hijij{0,LUi,), 

PGwK) = 16Tr^LO^GwlTTH,,,,{0,LJ^) , (21) 

where uji, is the angular frequency measured at the mo- 
ment of the gravitational waves generation. To obtain 
PGw('i^*) and the frequency today, one must account for 
the gravitational wave amplitude and frequency rescaling 
given by Eq. O 

Since equipartition between kinetic and magnetic en- 
ergy densities is maintained during Kolmogoroff turbu- 
lence, the total source for the gravitational waves (from 
the magnetic and kinetic turbulence) is simply two times 
the kinetic turbulence source. Using the assumptions 
made above, the source Hijij tensor is given by O, Ell , 



67r3/2 
erfc 



1/3^2/3 J ■ 



F2/3y^4/3 



(22) 



Here, erfc(x) is the complementary error function defined 
as erfc(a;) = 1 — erf(a;), where erf(x) = dye'}qi{—y^) is 
the error function, e — (ao/a^,)e and k = {ao/ai,)kQ 
are the physical energy dissipation rate and the physical 
wavenumber respectively. As can be expected, the inte- 
gral in Eq. (|22p is dominated by the large scale (k ~ fco) 
contribution so, for forward-cascade turbulence, the peak 
frequency is [13| Wmax.* — ^o^^o- It must be noted that 
the peak frequency is determined by the time character- 
istic of turbulence only in the case when the turbulence 
duration time is enough long {tt — few x Tp to insure the 
applicability of the Proudman argument 40] , which has 
been used to justify the use of the Kolmogoroff model. 
Otherwise, the peak frequency is determined by the char- 
acteristic scales of the pulse-like source [l^ . 

Taking into account the expansion of the Universe and 



using Eq. 



for the Hijij tensor, we find the gravi- 



tational wave amplitude as a function of the linear fre- 
quency measured today, / = (a*/ao)/* with fi, = uji,/2TT, 



1 /2 

X [TTUJi,HlHijij{Q,uu^)\ 



1/3 



where fn = is the Hubble frequency measured today, 
fn ^ 1.6 X 10-5Hz(g*/100)i/6(r*/100GeV) . 





FIG. 1: hc{f) for the EWPT with = 100, vt = 1/2, 
and vo = va with zero magnetic helicity. The top panel: 
T. = 80 GeV, (solid line), T, = 100 GeV (dash line), and 
= 150 GeV (dash-dot line) with /3 = lOOH*, while T* = 
100 GeV and (3 = 50H^, (solid line), = lOOH^, (dash line), 
and (3 = 150//* (dash-dot line) in the bottom panel. In both 
panels the bold solid line corresponds to the 1-year, 5a LISA 
design sensitivity curve [121 including confusion noise from 
white dwarf binaries, bold dash line [5311. 



(23) 



Eq. l[23l) allows us to predict the gravitational wave 
spectral properties: i) The low frequency (/ ^ /peak ) 
dependence is hc{f) oc f^^^, leading to Ogw(/) oc 
such a behavior is common for all causal sources [Tsl 
and it is true for any kind of the waves, including the 
sound waves generation by turbulence '54'|; ii) The peak 
position is determined by the time-duration of the source, 
and it is equal either to vaj/h (developed stationary 
source) flS*! or jfn (pulse like source) [ll|, [1^; iii) at 
higher frequencies, / ^ fa, the gravitational waves am- 
plitude is damped exponentially due to the exponential 
temporal decorrelation of fluctuations [l^ as opposed to 
the power law slow damping shape of the gravitational 
waves generated by the bubble collisions @, [a, [1, [l^ H^l ; 
iv) for the turbulence generated gravitational waves the 
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power law shape takes place in the vicinity of the peak, 
i.e. for /peak < / < Ri/^ff, hcif) (x dl. 

V. RESULTS AND DISCUSSION 

The results for hc{f) for the EWPT are shown in Fig. 
1. For a semi- analytical estimate it is straightforward to 
get the peak frequency of the amplitude of the gravita- 
tional waves emitted during direct-cascade to be (l^l 

= (^) (^) <^^' 

The peak frequency is shifted to the lower frequencies 
with T-t increasing due to Jh dependence on T^. Using 
the definition of the Hijij{0,uji,) tensor, Eq. (P^ . the 
gravitational wave signal reaches its maximal amplitude 
approximately at / = j~^VAfH, and then, 

According to Eqs. and (^5]) . the peak amphtude of 
the EW phase transition gravitational signal is order of 
5 X 10^21 for = 100 GeV, /3 = lOOH^, and = 100, 
with /jf(3c — 10^"^ Hz, which is an agreement with Fig. 1. 
For the stochastic gravitational waves the real LISA sen- 
sitivity will be lower. Even accounting for this, the grav- 
itational signal from the EW phase transition should be 
detectable, since it significantly exceeds the LISA noise 
around 2mIIz, if we adopt the model described here with 
vii = 1/2, and va — 0.3. Rewriting this result in terms of 
the gravitational wave spectral energy density parameter, 
we obtain Qcwif = /peak) ^ 2v%jHl00/g,p^ x IQ-^. 
Note i^Gwif = /peak) is the temperature independent 
for a given value of va and Vb, as it can be shown from 
the dimensional analysis and only slightly depends 
on 5*. 

The peak frequency for QCD phase transitions with 
n = 0.15 Gev, f3 = 6iJ*, and = 15 is 

G = X lO-VpcYk ^ 2 X 10-^ Hz. (26) 

which is order of six magnitudes lower than the LISA 
low-frequency sensitivity. 

We also define the efficiency of the gravitational wave 
production kqw as the ratio between the magnetic en- 
ergy density available from the phase transitions and the 
energy density converted into the gravitational waves, i.e. 



KQwif) = Pgw/pb- Since both total energy densities 
scale the same way with the expansion of the Universe, 
kgw is invariant (no damping of the magnetic field). It 
can be shown that kgw — 2(/3/i?*)~^w^wj^ <^ 1 [Ts*! (for 
EW phase transitions with the parameters mentioned 
above kqw — '^l^ips/ Prad)"^^^ — 10^®), so only the small 
fraction is transferred to the gravitational wave signal. 
However, the LISA sensitivity is ^Icwif — ImHz) ^ 
lO^^'^ and the magnetic field generated signal from the 
EW phase transition would be still detectable if even 1% 
of the radiation energy consists on the magnetic energy. 

We emphasize that although we present our results for 
Vh = 1/2, they can be easily extended for an arbitrary 
value of viy. For example, a different model of the EW 
MSSM phase transition gives Vb an order of magnitude 
smaller, see Ref. [11]. Another degree of the freedom 
is related to the energy scale of the phase transitions. 
In particular, for the EW phase transition the energy 
scale is approximately equal to the Higgs mass, which 
ranges from about 110 GeV to 127 GeV [i^, however, 
such freedom of the EW phase transition energy scale 
just slightly affects our final results. 



VI. CONCLUSIONS 

The model presented above might be applied for dif- 
ferent mechanisms [1^ of the magnetic field genera- 
tion leading to the substantial magnetic energy presence 
during phase transitions, but emphasize that the ours is 
based on the derivation of the magnetic field amplitude 
from the fundamental EW MSSM or QCD Lagrangians. 

Summarizing, using MHD turbulence model with no 
helicity we find that the gravitational wave produced dur- 
ing the EW phase transitions is most likely detectable by 
LISA, while that produced by the QCD phase transitions 
will not be detectable. 
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